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^SJ Abstract 

O We consider the Fast Diffusion Equation Ut = Au'^ posed in a bounded smooth domain f2 C R'* 

with homogeneous Dirichlet conditions; the exponent range is iris = (d — 2) + /(d + 2) < m < 1. 
It is known that bounded positive solutions u{t, x) of such problem extinguish in a finite time T, 
and also that such solutions approach a separate variable solution u(t,x) ~ (T — t)^'''^~'"'S'(a;), as 
t — >■ T~ . Here we are interested in describing the behaviour of the solutions near the extinction 
I— I time. We first show that the convergence u{t,x) (T — to ^(a;) takes place uniformly in 

Qi^ the relative error norm. Then, we study the question of rates of convergence of the rescaled flow. 

For m close to 1 we get such rates by means of entropy methods and weighted Poincare inequalities. 
• The analysis of the latter point makes an essential use of fine properties of the associated stationary 

elliptic problem — A5"" = in the limit m — >■ 1, and such a study has an independent interest. 
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1 Introduction 



We are interested in describing the behaviour of nonnegative solutions of the Fast Diffusion Equation 
(FDE) near the extinction time. More precisely, we consider the following initial and boundary value 
problem 

{M^ = A(u'") in(0,+oo)xf7 
u{t, x) — for T > and x € dfl , 

posed in a bounded connected domain C M'^ with a regular boundary of class C^'", a > 0. The fast 
diffusion range is < m < 1, but the theory developed below needs the further restriction < m < 1, 
where the lower end is the exponent nis — [d — 2)+/{d + 2) (inverse Sobolev exponent). We assume 
that the initial data uq is a bounded and nonnegative function. It is well-known that the above problem 
possesses a unique weak solution u > that is defined and positive for some time interval < r < T 
and vanishes at T = T{m,d,UQ) > 0, which is called the (finite) extinction time, cf. [71 [TU [511 US]- 
Note that the conditions on the initial data can be relaxed into uq g LP(ri) for some p > 1 and p > pc 
where Pc = max{l,(i(l — Tn)/2} in view of the LP-L°° smoothing effect, see [Hll^ . 

Rescaled equations and previous results. To study the asymptotic behaviour it is convenient to 
transform the above problem by the known method of rescaling and time transformation. Thus, we put 

(1.2) u{t,x)= f^^^)' ^ v{t,x), t = T\og^ ^ 



T ' " °\T-T 



In this way. Problem (1.1 1 is mapped into the equivalent "rescaled problem": 



vt = A(t;™) + in (0, +oo) x ft, 

(1 — mji 

(1.3) { v{0,x)^uo{x) infl, 

v{t, x) =0 for t > and x e dil. 

The transformation can also be expressed as 

(1.4) v{t,x) u(t -Te-*/^,xY 

and the time interval < t < T becomes < i < oo, so behaviour near extincion for the original flow 
becomes behaviour as i — > cxd in the rescaled flow, which is more convenient to analyze. Thus, in a 
celebrated paper, Berryman and Holland [7], 1980, reduced the study of the behaviour near T of the 



solutions of Problem ( 1.1 1 to the study of the possible stabilization of the solutions of the transformed 



evolution problem (1.3 1. They showed that the solutions of the latter problem stabilize towards the 



solutions of the associated stationary problem 

(1.5) - A(S"") = J^^^^ S in n, S{x) =0 for x e dn , 

where < m < 1 and il are as before, and 5 > in fJ. Using the new variable V = S"™ and putting 
p = 1/m > 1 and c = 1/((1 — m)T) the latter problem can be written in the more popular semilinear 
elliptic form 

(1.6) -AV = cVP in n, V = on dn 

Note that our restriction m > is the exact condition that makes the last problem subcritical, 
P<Ps:^ (d + 2)/(d-2). 
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It is precisely proved in [T that the rescaled orbit of a solution v{t) — v{-,t) converges in M^Q^'^(ri) 
along subsequences to one or several stationary states S (Remark: the elliptic problem can have multiple 
solutions depending on the geometry of 51). In the language of dynamical systems, the omega-limit of 



V is included in the set of positive classical solutions to the stationary problem (1.5 1. 

A main issue that remained open after [7] was to understand wether the rescaled solution v converges 
to a unique stationary profile for alH — )■ oo, even when the set of stationary solutions contains more than 
one function. The question of uniqueness of the asymptotic profile has been solved by Feireisl-Simondon 
in |25) . We rewrite their main result in our notations: 

Theorem 1.1 ([25]) Let fl C M'^, d > 1 be a bounded domain of class C^'", a > 0. Let v G 



L°°((0, oo) X il) be a bounded weak solution of (|1.3|, then v is continuous for t > and there ex- 



ist a classical solution S of the stationary problem (1.5) such that 



v(t) — > 5 in C(Q) as t oo . 

We recall that Theorem 3.1 of [25j is a bit more general, indeed we specialize here to the case f{u) 
—uj (1 — m). 



We remark that every solution S — Sm.rix) to the elliptic problem (1.5) produces a separable solution 
U of the original FDE of the form 



(1.7) U{t,x)^S{x) 



T-T 

T 



which corresponds to the initial datum l/({0,x) = S{x). Indeed, it is not a solution but a family of 
solutions since we can fix T > at will: we will write U t for definiteness when needed. 

In the present paper we present improvements on these results in several directions: 

(i) We prove that the stabilization process takes place with Convergence in Relative Error. This topic 
occupies Section [2] and the main result is Theorem 2.1 



(ii) In Sections [S] and [5] we prove convergence with rates to the stationary state using so-called entropy 
methods. The use of the term entropy deserves an explanation: we introduce a suitable Lyapunov func- 
tional, namely a weighted L^-norm of the quantity v(t) — v{t), that decreases in time along the nonlinear 
flow, and moreover, its dissipation in time is carefully controlled. We have decided to call entropy such 
functional since it has properties similar to the entropy functionals that have been extensively used to 
study the Cauchy problem, in particular the one used in [3191 [10]; we do not claim a physical meaning 
for the entropy that we use here. 

The first step is contained in Section |3j where we obtain convergence whenever a certain Weighted 
Poincare Inequality holds with a sufficiently large constant, and precise decay is shown in that case, cf. 
Theorems 13.31 and 13.51 



(iii) Then, in Section [5] such assumption is shown to hold for the solutions of our problem in a restricted 



exponent range m# < m < 1, and we obtain the concrete asymptotic results. Theorems 5.6 and 5.17 
The study relies heavily on the analysis of the associated semilinear equation ( |1.6| ) in the limit p ^ 1 
which has in our view an independent interest and is developed in a previous Section |4[ The main result 



is Theorem 4.1 In it we show that, if is a solution to equation (1.6 1 with homogeneous Dirichlet 
boundary datum and we choose c = Ap > in such a way that ||V^||p+i = 1, then Vp/^i — > 1 uniformly 
in as p — > 1, where $1 is the ground state eigenfunction of the Dirichlet Laplacian on fi, with the 



normalization ||$i||2 — 1- Moreover we have that Ap — ^ Ai as p — > 1 and, finally, Proposition 5.4 also 
shows that 1/(1 — m)T — )■ Ai as m — >■ 1. 

(iv) The entropy method applies also for the Porous Medium Equation, that is when m > 1, and it 
allows to find the rate of convergence, thus recovering the sharp result of Aronson and Peletier [3], by 
different methods. We devote Section |2.4| to present the the slow diffusion case m > 1. It is worth 
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mentioning that the we also obtain, as an intermediate resuh, a faster convergence for the entropy 
functional used in Section [3j see Theorem 3.4 which is new. 



Notations. Before proceeding with the statement and proofs of the results, let us recall some nota- 
tions. S will denote the stationary solution of Problem |1.5| indicated by Theorem For x G ft we 
write d{x) — dist(a;, dil) to indicate the distance to the boundary, its properties will de described below 
when needed. Ai is the first eigenvalue of the Laplacian operator in the domain with zero bound- 
ary conditions and positive eigenfunction $i. 52 is the optimal constant in the Sobolev embedding 
Wa'^ifl) 2* = 2d/{d - 2), d > 3. For d > 3 we put = 2* - 1 ^ {d + 2)/{d - 2) and 

iris — {d — 2)/((i + 2) — \/ps- By || • ||p we denote the standard L^'(r2)-norm, 1 < p < oo, other norms 
will be carefully denoted. 



2 Convergence in relative error 

We show that the quotient v/S converges to 1 uniformly in the whole of f2, up to the boundary; 

Theorem 2.1 Let u he the solution to Problem \l.l\ and let T — T{m, d, uq) be its extinction time. Let 
S{x) be the positive classical solution to the elliptic problem (1.5) predicted in Theorem \l.l\ Then, 



(2.1) 



lim 

T^T- 



u(t, •) 



Uir,-) 



lim 

i— >-oo 



Si-) 



where U is the separable solution (1.7) and v is rescaled solution defined in (1.4 1. 

This type of convergence is what we call uniform relative-error convergence (REC for short), and it 
is our first main contribution to the subject of fine asymptotics near extinction. We can rephrase the 
result in terms of the rescaled solution v as 



(2.2) 



lim 



S"-{t,-) 



- 1 



= 0. 



L°°(0) 



and this form will be practical for the calculations. Since V — 5™ is a function with linear growth near 
the boundary, what we say is that 



(2.3) 



[v"'it,x)- S"'ix)]/dix)^0 



uniformly in x G 17 as t — >■ cxd. 

As for related results, DiBenedetto, Kwong and Vespri proved in [21], the Global Harnack Principle 
that we recall with our notations as follows, since we will be using it throughout this section: 



Theorem 2.2 ([24]) Let v be the solution to the problem (1.3) with rus < m < 1. Then, for any e > 
there exist positive constants Ca,m , Ci,m > depending on d, m, ||uo||m+i, IIVm^Hj, dfl and e, such 
that for any t > e and for any x G Q 

(2.4) Co,™ dix) ^ < v{t, x) < Ci.m d{x) ^ , 

Moreover for every k > 0, t > e and for any a; G 



(2.5) 



Mt,-)\\c-(n) = max sup \D°'v{t,x)\ < ' ' <a;)- 
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Note that the constants Co,m , Ci,r?x > may degenerate as to — > 1 or m — > mg. Estimate (2.4) above 



immediately imphes the foUowing estimate for the sohition to the FDE in the original variables: 
(2.6) Co,™d(a;)i/"(r-r)i/(i-™) < «(r,x) < d(a;)i/™(r - r)i/(i-'"). 



This result is weaker than our Theorem |2.1| in the sense that our relative error convergence result is not 
a consequence of the above estimates, neither in the original nor in the rescaled variables: inequality 



( 2.6 1 only implies that the quotient v(t)/S is bounded and bounded away from zero up to the boundary, 
but it does not prove that it converges to 1 as t — > oo. As far as we know, only the papers [71 [251 [2^1^ 
contribute to the subject of the asymptotic of the Dirichlet problem for the FDE on bounded domains. 

On the other hand, there is an extensive literature on stabilization of solutions of evolution equations in 
different norms, mainly in or C" spaces. Let us comment on some results on the topic of convergence 
in relative error which are not so usual. Uniform convergence in relative error was first proved for the 
Fast Diffusion Equation by one of the authors "41] in the following setting: solutions are nonnegative 
and the equation is posed in the whole space R'^ with exponents < to < 1, where m,, = {d — 2)^/d. 
The case of all to < 1 is treated in [51 [TTJ [TUI where sharp rates of convergence in relative error are 
obtained for FDE posed in IR'* with a new entropy method, that does not apply to the bounded domain 
case. 

In the present setting of homogeneous Dirichlet data in a bounded domain, the result is true for the 



Heat Equation (see a brief account in Section 3.1). It is also true for the Porous Medium Equation, 



i.e., our problem for to > 1, and the result follows from analyzing the asymptotic result of Aronson and 
Peletier [3] ; see the survey paper [42] . We recover the known results also in the PME case, which turns 
out to be simpler, see Section |2.4| . Hence, our present interest in the Fast Diffusion case, where the 
presence of extinction makes the boundary argument more difficult, since the usual super-sub solution 
method does not work. 

In the next subsections we proceed with the proof of Theorem |2.1| 



2.1 The relative error function and its equation 



From now on we will consider the evolution problem in its rescaled form (1.3). Theorem |1.1| proves 
that the w-limit of a solution is contained in the set of classical solutions to the Elliptic Problem ( |1.5[ ), 
and the convergence takes place in the uniform norm, and the solution v selects a unique profile S to 
converge to. Let us fix it: once we consider v{0,t) = uq then we know by Theorem 1 1.1 [ that v{t) — ^ S. 

Now we introduce the Relative Error Function (REF) 



(2.7) 



V" 



, + i) = y(0 + i), 



and 



• Tiie Parabolic Equation of the REF and the regularity of its solutions. Using the equations satisfied 
by V and v, is easy to show that cf) satisfies the following parabolic equation 



(2.8) 

where F is given by 
(2.9) 



1 



(1 



*t = V'-^Acj) + 2^ • V0 + F(0) 



F(</.)-c[(l + 0)i/™-(l + , 



Estimates (2.4) on v and S, which is a stationary solution, imply that 



< 1- 



— C2,m < < C^_m — 



G 



0,m 



- 1 
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which proves that (j) is bounded uniformly in (t, x), for t > to > 0; notice that to can be chosen arbitrarily 
small, but this affects the value of the positive constants Co.™ and Ci,„i. Moreover, we notice that 

l + (P^v"'/V>0 

in the interior of il. Since (j) is also bounded in the interior of Q, we conclude that the parabolic equation 



(2.8) is neither degenerate nor singular in the interior of fl. It follows from standard quasilinear theory 
(cf. [32]) that the solution cj) of such a parabolic equation is Holder continuous in any inner region 
r^/ C (the fact that (j> is Holder continuous could also be proved by observing that (p + 1 = 
and by recalling that both v (see e.g. [Mj) and S (see e.g. [23] or [7]) are at least Holder continuous 
and positive in the interior of fl). 

• Convergence of the REF in an interior region of fl. Under the running assumptions, we know by 
Theorem O that 

sup \v{t) — 5| — > as t ^ CO , 
n 

but this is not sufficient to prove the convergence of the quotient v"'' / S"^ to 1 in the whole il, since at 
the boundary there is the problem caused by the fact that both v and V are zero and that the parabolic 



equation (2.8) can degenerate at the boundary. However, such a problem is avoided in any interior 



region where both v and S are strictly positive. We define such interior region as 

^i,s ^ {x eQ : dist(a;, dQ) > 6} , 

with distance from the boundary (5 > which will be chosen later small enough; we can thus say that 
in any interior region 51/ ,5 C we have 

ll0(*)llL-(o^,i) = sup |(?;'(t, •)! ^ 0, as t-^oo 

recalling that </> = t;™/S"" — 1. We can sum up what we have proved so far in the following Lemma 
on Inner Convergence: 



Lemma 2.3 Let w be a solution to the rescaled problem (1.3), and let (p be the corresponding relative 



error function defined by (2.7). Then for every £ > and 5 > Q there exists to, £,5 > 0, such that for 



every t > to.e.a cind for every x G ^1^5 we have 

(2.10) m,x)\<e. 



2.2 Distance to the boundary and barriers 

To get the proof of the convergence theorem we still have to show that uniform convergence of (j) takes 
place up to the boundary. To this end we will use a barrier argument, based on the following lemmas. 
We remark here once and for all, that the barriers are independent of the particular choice of the 
stationary solution S. 

First, we collect some properties of the function "distance to the boundary". It is defined as usual: 

d{x) — dist{x, dil) ^ inf |a; — ?/| 

where | • | is the Euclidean norm of M"*. 

Lemma 2.4 (Properties of the distance to the boundary) Let fl C M.'^ be a bounded domain 
with boundary dfl of class . Let for ^ > 

Q.^^{x(^Q. : d{x) <£}^VL\ Oj^ 
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be the open strip of width ^ near the boundary. Then, 

(a) there exists a constant > such that for every x G f2^(, , there is a unique h{x) G 9f2 which realizes 
the distance: 

d{x) = \x ~ h{x)\. 

Moreover, d{x) G C^(f2^j,) and for all r G [0,^o) the function : d{Q.r) n 57 — J- dVl defined by 
Hr{x) — h{x) is a homeomorphism. 

(b) Function d{x) is Lipschitz with constant 1. i.e. 

\d{x)-diy)\ < \x-y\. 

Moreover, 

< c < |Vd(x)| < 1, for any x G fl^^ 
and there exist a constant K > such that: 

(2.11) - K < Ad{x) < K, for any x G Q^g 

We refer to for the proof of this lemma. Part (a) is due to Serrin. 

We need a second technical result about estimates for the gradient of the function V = S*™ near the 
boundary. 5 is a positive classical solution to the elliptic problem (1.5) 

Lemma 2.5 For > small enough there exists /?o > such that 

vv{x) ■ vd{x) > /3o > , yxen^„. 



Proof. As explained in Section 6.2 or as a consequence of the bounds (2.4), the following estimate 
holds 

(2.12) C^,nd{x) < V{x) < Ci"„,d(x), for any x e Q . 

Moreover, we know that S* is a positive classical solution of the elliptic Dirichlet problem |1.5| since 
TTT-s < m < 1. Thanks to Lemma |2.4| we can conclude that both V and the distance function d are 
functions of class in a suitable neighborhood of the boundary fi^g . The above estimates imply 

^rn d{x)-d{xa) V{x) - V{xo) d{x) - djxo) 

- \x-xo\ - \x-xo\ \x-xo\ ' 

for any a; G f?^ and any xq G dfl, since d{xo) — V{xq) — 0, if xq G dfl. This implies that 

< C^M^o) < d,V{xo) < CZndA^o) 

and that djV{x) and djd{x) are both nonnegative, so that WV{x) ■ \/d{x) > 0. Moreover it implies 
that: 

3 3 3 

and finally: 

(2.13) < cc^-^ < Co";„jvd(xo)|2 < vv^(xo) • vd(xo) < cr;„|vd(xo)|2 < CT,^ 

for any xq G dVl, since we know from lemma ( |2.4[ ) that 

< c < \Vd{x)\^ = ^{djd{x)f < 1, for any x G . 
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By continuity of Vd(x) and since |A(i(a;)| < K, we can extend the estimates (2.13) from xq E dil, to 
a small neighborhood of the boundary, say x € fJjo, eventually by putting a smaller lower constant 
< /3o < cC^,„, which can eventually depend on if > and • □ 

We have obtained uniform estimates on a neighborhood of the boundary fl^^ . These lemmas are needed 
for a key ingredient in the proof of the relative error convergence theorems, which is the construction 
of barriers as super-solutions. 

Lemma 2.6 We can choose positive constants A, B , C so that for every to > the function 
(2.14) ^{t,x)^C -Bd{x)- A{t-tQ) 



is a super- solution to Equation (2.81 (the one satisfied by the REF (p) on a parabolic region 

= _i = {(<, x) e (<o, oo) X f7 : $(t, x) > -1} , 



and moreover T,^ C (to,oo) x il^j^. Super-solution means that for any {t,x) G S]$ we have 
- (1 + $t > T/1-i/™a$ + ■ V$ + 

with 



1 



1 — rn 



(l + $)^/"-(l + $) 



> 



and = S*™, S* 6em^ a positive classical solution to the elliptic problem (1.5). The constants A,B,C 
depend only on m, d, the upper bound for and the geometry of the border through Co^m cind Ci^m- 

Remark. The construction of the barrier is quite technical, so we stress that when is small enough, 
then a sufficient condition on the parameters is 



(2.15) 



1 — ~ 777 

1 + C+ A] ^i<{l3By 



We considered the barrier on the parabolic region since in that region the quantity $ + 1 > 0, 

but in what follows we will only need the smaller region S^.e , for small e > 0. 

Proof. We recall that the equation satisfied by the REF is 



(2.16) 



m 



vl/m-l 



yi-i/™A^ + 2_.V0 + F(0) 



We will prove that function (2.14) is a super-solution for equation (2.16) on the parabolic region I]$ if 
we find constants A, B and C such that 

(2.17) + $t > (/) > (//) > 1/1-1/™ A$ ^ ■ V$ + F($) 

where the expressions (I) and (II) are estimates from below and above respectively for the left and the 
right terms, independently of {t,x) € S$. 

(I): This estimate is simple since $t = — A for any (t, x) £ E$. Hence: 



-(l + $)i/™-icO, = --(l + $) 



l/m-l 



{!)■ 



(II): This estimate is more involved. First, we rewrite the right-hand side of (2.171 in a more convenient 
form: 



F1-i/™A$ + 2 



V$ + F($) 



-B 



[V Ad-^2W -Vd] +F($) 



< 



B fl + $ 

[VAd + 2\/V-\/d]^ ^ 



l/m 



1 — m 
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since 



A$ = -BA{d{x)), V$ = -BV{d{x)), 



1 



< 



1 — rn 



($ + 1) 



1 — m L 
Moreover, we have that 

(2.18) [V Ad + 2VF • Vd] > -XCJ';,,,^ + 2/3o = /? > 

on a region for any x G 51^^ with 



(2.19) 



^1 = min <^ ^0, 



2/3o 



with if, /5 Ci.m > as in the previous two lemmas. Indeed, we have 

V{x)Ad{x) > -KV{x) > -KC'i^rn d{x) > -KC'i^rnii 



for any x G fi^^, as a consequence of estimate (2.11) and estimates (2.12). Moreover, 

VV{x) ■ Vd{x) > /3o > 



for any x € 51^^ , as proved in Lemma 2.5 FinaUy 



y-v^Ac + 2^ . va> + F($) < + il±^ 

yl/m \ ^ - 1-7] 



1 /t?j. 



in fi^j with ^1 > as in (2.19). With these estimates, we conclude that $ is a super-solution if the 
following condition holds: 



A 

m 



(1 + $)i/™-i = (/) > (//) = -P+ (1 + $). 

^ ' \ J - \ J 1 - m 



We can rewrite it in the form 
(2.20) 



A 1 + $ 

- + 1 

m 1—771 



(i+$)i/™-i < 4^ 

^ ^ — ^l/m 



with (t, x) e E$ and > as in (2.19) 



We have thus proved that <I> is a super-solution for any {t,x) G with x € Sl^^ such that (2.20) 
holds. Since $(i,x) < C in the region under consideration, it suffices to choose A, B and C so that 



(2.21) 

to ensure that (//) < (/). □ 



A l + C 



m 1 — m 



(i + c)i/"-i < 



2.3 Proof of Theorem [2J] 

It is based on the previous Lemmas and Theorem |1.1| of |25) . 

(I) We have to show that given e > there exists a time T(e) > such that for any t > T{e) and for 



any x £ Q we have that |0(t, a::)| < e. By Lemma 2.3 we know that for every e > and S > there 
exists io,e,<5 > 0, such that for every t > io.e.a and for every x in interior region ni,s we have 



\(f>{t,x)\ < e. 
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It remains to show that uniform convergence takes place also up to the boundary. This will be a 
consequence of comparison with the barrier function of Lemma |2.6[ given by 



x)^C - B d{x) - A{t - to). 

A, B, C of the barrier $ are suitable positive constants which are chosen as in Lemma [276} while tQ > Q 
is a free parameter which will be adjusted later. 

(II) Comparison of (p with $ takes place in a neighborhood of the parabolic boundary of the form 
Q* = (toiT) X fls for 6 small, and tQ,ti to be determined. The parabolic border of this region is 
formed by three pieces: the initial section at t = to, the inner parabolic boundary, and the outer lateral 
boundary. In order to compare (p &nd <I> we have to check their values on the above three pieces of 
parabolic boundary. 

(a) We compare the values of cf) and $ at the initial section t — to- We want that 

Hto, x) < $(io, x) = C - Bd{x) 
for all X E ils- This is possible because of the uniform boundedness of (j) 

C2,m = c{Co,m — Ci^m) < 4>{toTx) < C^^m = c(Ci^m — CQ,m) 

for all a; e as a consequence of bounds (|2.4[). Now we simply choose C sufficiently large previous to 



the choice of A and B that have to satisfy (2.15) 



(b) Comparison on the inner parabolic boundary: This piece of the boundary is given by the points 
{x,t) such that d{x) = 6 and t e {to,T). On this set we want (l){t,x) < ^{t,x). Let us fix e > and 



< i5 < ^1 where > is given in Lemma |2.6| By the uniform inner convergence (cf. Lemma 2.3 ) 
we know that there exists a t*{e,S) > such that (l){t,x) < £ on ilj^s if t > t* . The desired comparison 
holds if 

(2.22) e<C-B5-A{t-to). 



Since C cannot be small this implies restriction on B that has to be compatible with (2.15). This 
happens if S is small enough. Once B and C are chosen, it suffices to take A{t — to) small. 

(c) We still have to check the comparison at the outer lateral boundary, [to,T] x dH., where we only 
know that (j) = v"^ / S"^ — 1 is bounded. But we can use an approximation trick using the solutions u,, 
of problems posed in the domain fl^ which is smaller than fl. We know that u,, w as ?y — > 0. 

Using Urj instead of u and Qri instead of allows to say that 0^ = u.^/lA — 1 = on the new outer 
boundary, hence (pr, < 0. Then we have: 

(pr, = < B5 + e < C ~ A{t ~ t*) = <!> 

thus (/),,< $ also on the outer boundary dflrj- 

Parabolic comparison allows then to say that < $ in the region Q* for t > to > t* such that 
t — to < {C — BS — e)l A. Pass to the limit in 77 ^ to get < $ in Q*. In this way the following 
improvement of convergence near the boundary after some time delay given by 

C-B8-e 

he,& = 1 , 



which is the maximum that (2.22) allows 



Steps (I) and (II) can be summarized in the following 
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Figure 1 : Idea of the behaviour of the harriers: 
y-axis : values of^{t,x) 

X-axis: values of d{x) = d{x,dVl), i.e. the distance from the boundary. 

Sj ; the points where '^{t,x) — Ei, i.e. the points of the boundary , z = 1,2,3. 

Si : different values of e (decreasing with i = 1, 2, 3 J give different barriers 

decreasing with e as the arrow (i) indicates. 
^1 and 5 are as in Lemmas \2.6\ and \2.7\ 



Lemma 2.7 Under the above conditions we have for t = to + hi^^s 
(2.23) 4>{t,x)< 



e + BS, 



for any x d fl, such that d{x,d^) > S, 
for any x € fl, such that d{x,dfl) < S. 



provided that to > t* 



(III) The proof of Theorem 2.1 in the version of formula (2.2) follows now by fixing e > 0, finding a 
barrier with constants A, B, C and then taking S < e/B. If t > t*{e, S) + h,, s, then 

(t){t,x) < 2e 



everywhere in Q. □ 



2.4 The Porous Medium case 

In this section we consider for the sake of comparison the same Dirichlet problem for the Porous Medium 
Equation 

u^ = A(u") in(0,+oo)xf7 
u(0, x) = uo{x) in n 
u{t, x) = for r > and x S d^l , 

now for m > 1. In this case there is no extinction in finite time, which changes the previous analysis 
and makes it much simpler. Let us give some details. By means of the rescaling 

u(t, x) — u{t, x){1 + r) ™-i and 1 + r = e* , 
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the problem is mapped into the equivalent "rescaled problem" 

V 



m — 1 



vt = A(t;'" 
v{0, x) — Uq{x) 
v{t,x) = 

The transformation can also be expressed as 

v{t, x) = e™-! u (e* — 1, x) 



in (0, +oo) X ri, 
in ri, 

for t > and x E dfl. 



with 



t - log(l + r) 



and the time interval < r < cx) remains < t < oo, in particular we preserve the initial datum. 
Notice that the rescaled problem when m > 1 is formally the same as for the Fast Diffusion case, if 
one just considers c = — 1) > 0. In this context it has been proved that solutions converge to a 

stationary state S, which is the unique solution of the elliptic equation —AS*™ = cS", see for example 
[21 142]. Solutions by separation of variables are given by Z^(t, x) — S{x){l + r)^^/^™^^^ . The optimal 
rate of convergence in relative error as r —> cxd 







W(r,.) 


L°°(n) 



u(t,-)(1 + t)™-i 



Si- 



< 



K 



L°°(n) 



1 



has been first obtained by Aronson and Pcletier in [3], Theorem 3, for smooth nonnegative initial data. 
The rate can be easily shown to be optimal because of the special family of global solutions 



Uk{T,x) 



Six) 



for any fc > . 



Indeed, the proof of [3] strongly uses comparison with the stationary state S through the special solutions 
L(k- In the fast diffusion case comparison with smaller solutions does not help, since they extinguish 
earlier. This is just a example of the extra difficulties in proving convergence in relative error when 
TO < 1. In rescaled variables the result of 13 reads 



(2.24) 



\v{t,x) - S{x) \ < KS{x)e 



for all a; e and t ^ 1 , 



Other interesting approaches to convergence results, together with extensions to a larger class of initial 
data and solutions, can be found in [42]. Our entropy method applies to the case to > 1 as we shall 
briefly discuss in the next sections, and allow us to recover the result of [42] and the optimal rates of 



convergence in relative error of J3 , cf. Theorem 5.8 



3 Stabilization with convergence rates 

We recall here the setup and the notations. We consider the rescaled equation 

= A(w'") + cw, 



(3.1) 



1 



c = 



r(i - to) 



> if TO < 1, and 



1 



TO — 1 



> if TO > 1 



posed in a bounded connected domain 17 C M"^ with sufficiently smooth boundary. We mainly deal with 
the so-called fast diffusion exponents and then we assume to^ < to < 1. Almost all the calculations will 
hold for any to > to^, including the case to > 1; we will emphasize the differences when they will occur. 
The linear case to = 1 is well known and will be briefly recalled in the next subsection as a motivation 
of our techniques. 
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We now introduce the quotient w — v/S, which converges uniformly to 1 on fl, as a consequence of 
Theorem |2JJ We then have 



(3.2) 



Swt = A(5'"w") 



where 5' is a stationary solution so that A5™ + c5 = in il with = on dft, as precisely indicated 
in Theorem [lT| We will also write V = 5", which satisfies AV + cVP = in fi, = on with 
exponent p = l/m > 1. 

We propose to perform the calculation on the asymptotic decay in terms of^^w — l = u/S'— 1, 
which is the relative error of the solution. Notice that it is different from the relative error (j) used in 
the previous section, formula (2.71, which was defined a.s (p = — 1, hence (p = {1 + 0)™ — 1. We have 

1 



(3.3) 0t = -A(^'"(l + 0)'")+c(l + 

Using the identity 

A(5"(l + ffD) = V • [5'"V(1 + 0)™] + V(5'") • V(l 
and the equation for S*, this can be further written as 



AS"" (1 + ey 



1 



V • (S""V(1 



■v{i + er + cf{e), 



" s ^ y ■ , > g 

where 

(3.4) fie) := {i + e)-{i + ey' = (i - m) + ^ q^^s^ 

for small 9. One more calculation gives the following form for the equation 

(3.5) S"'+^er = V • (S'2™V(1 + 61)™) + c f{e) 
that we will use below. In terms of V ~ 5™ and p — I /m we have 

(3.6) VP+^9t = V • (V'^Vil + e)"") + c yP+i f{9) . 



3.1 Weighted inequalities for the hnear heat flow 

In search for inspiration on how to proceed further, we compare the situation with the standard way of 
treating linear equation = Au, which has striking formal similarities even if there is no extinction 
in finite time. After the suitable linear rescaling, which takes in the linear case the form v{x,t) = 



e^^*w(x,t), we arrive at the same equation (3.1) with m = 1 and c = Ai. The role of the stationary 
solution S is now played by the first eigenfuction <I>i > of the Laplacian with zero boundary conditions 
(or any of its nonnegative multiples). The equation for 9 — v/^i — 1 is: 

(3.7) 9t = A9 + 2V<i>i-V9 = <S>-^V-{<PlV9), i.e., 'i>l 9t ^ V ■ {<PIV 9) , 



to be compared with (3.5). We recall next the known way to proceed with the asymptotic analysis of 
this equation via an Intrinsic Poincare inequality. Indeed, let us consider the self-adjoint operator —A 
with Dirichlet boundary conditions on dfl. This can be defined by a general procedure by defining it 
to be the unique self-adjoint operator associated with the closure of the quadratic form |Vupda;/2, 
initially defined for u G C^{il). Such operator has purely discrete spectrum, and we denote by > 0, 
J = 1, 2, . . . its eigenvalues, arranged in nondecreasing order, and by $j the corresponding L^-normalized 
eigenf unctions. 
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The spectral representation for the corresponding heat semigroup Ur = Au shows that, letting uq be 
the initial datum and Cj = J^^ uq^j da;, one has 

oo 

u{x,t) = Cje'^^^^j{x) 



so that 



C2 $ 



2g-(A2-Ai)t^ 



cie^^i*$i t-i-+oo ci $1 

In other words, the solution u{t) is close to the explicit solution Ui{x,t) = cie^'^i*<i>i and the relative 
error w between such solutions, defined above, decays exponentially in time with a rate A2 — Ai. Notice 
in addition that the spatial factor ^2/^1 is bounded. 

To prepare the way to recovering a result of this kind in the nonlinear setting, where no spectral 
representation is available, we reformulate the above property as follows. Starting from equation (3.7), 
it is then natural to investigate the behaviour of 9 by working in the weighted space L^(<i>^ dx). In fact 
first we observe that the weighted mean is preserved: 



_d 
dt 



/ 9^ldx^ [ V • ($?V6i) dx = 0. 
Jn Jn 



Then we notice that: 



d 

~di 



[ 9^<Pldx = 2 [ 6iV-($?V6') dx = -2 [ \V9f<^jdx. 
Jn Jn Jn 



By the the above conservation of weighted mean we can and shall assume that 0$^ = 0, where 



5*1 



In^ldx 



It is then clear that in order to get a decay rate for E[9] ~ J^^ 9^^\ dx it suffices to prove the following 
intrinsic Poincare inequality: 

Proposition 3.1 Let f € WQ''^{fl) and g = Then the following inequality holds 

(3.8) (A2-A1) / |5-5*j'$?dx< / \Vg\^<S>ldx. 

Jn Jn 

Although this inequality is well-known, we provide in Appendix [6]2 a short proof for the reader's 
convenience, and we recall there some sharp upper and lower bounds on the spectral gap A2 — Ai. 



3.2 Energy analysis of the nonlinear flow 

Inspired by the preceding linear analysis and after carefully choosing among the different options to at- 
tack the nonlinearities of our evolution process, we are going to prove a certain type of entropy/entropy- 
production inequalities. We define the suitable entropy functional to be 

(3.9) 8[9{t)] = \j^ \9{t) - m\^Sl%T dx , 



where 

m 



S,,9{t,x)Sl+^dx 
Jn '-'^'^ '^•^ 
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and S ~ Sc,m is the chosen solution to the eUiptic problem 
(3.10) 



S >0 iiifl 
S = ondn. 



c > and rus < m < 1, and the relative error satisfies the equation 
1 



(3.11) 



5*" 



:^V- (5"'"V(l + 0)") + c/(0), with /(0) = (1 + 0) - (1 + 1 



We do not specify boundary conditions, but we know that 9 is contin uous up to the boundary dQ and 
the convergence in relative error valid in C(f2), proved in Theorem 2.1 indicates that the boundary 



conditions stabilize to as t — >■ oo. It is not restrictive to assume that \9\ < e on 957, or by the 
maximum principle in the whole fi, for arbitrarily small e > 0. The price that we have to pay is just a 
time shift. We now prove the following 



Proposition 3.2 Let c > be as in (3.11, and nig < m < 1. Let 9 be a global smooth solution to 
equation (3.11) and let e{t) := \\9{t, •)||oo 0. Then, the following inequality holds 

(3.12) 



-:^£[9it)] >m[l + £(<)]"~^ / \V9{t,x)\' S'"'dx-2c[l-m + eit)]£[9{t)] 

for all times t > to where ta is such that e{t) ~ \\9{t^ OIloo < 1 for all t > tg. When m = 1 we recover 
the standard L'^ -weighted estimates that hold in the linear case, see Subsection 3.1 Moreover, when 
m > \, the property e(t) holds and, moreover, for t sufficiently large: 



(3.13) 



_d 
'dt 



£[9{t)] > m[l ~ £(<)]"-i / |V6lp 5*2" da; + 2c[m - 1 - e{t)]£[9it)] > . 

Jn 



Proof. Uniform convergence in relative error holds as t -> oo by Theorem 2.1 that means \\9{t, ■)\\oo ^ 
0. We set now d/i — S^^^ dx and we will write S — Sc,m throughout the proof, since no confusion will 
arise. First we notice that 



(3.14) 

and also that 
(3.15) 



[9{t) - 9it)] [dt9it)] dfi = dt9{t) / [9it) - 9{t)] d/i = 



0< / \9-9{t)\d^i^ / [9^ ~9itf]di^. 
Jn Jn 



We next differentiate £[9{t)] along the flow 
d 



.£[9it)]^ / [9{t)~9it)] [dt9it)]dti+ / [9it) - 9{t)] [^Mt)]d^l 

= - / V [9[t) - e{t)\ ■ [52™V(1 + 0)™] dx + c ( [9{t) - 9{t)\ f{9) dy. 
Jn Jn 

= -m /" (1 + 9)"'-^ |V6i|2 dx + c [ [9{t) - 9{t)] f{9) d^i = (/) + (//) . 
Jn Jn 

In order to estimate (//) we notice that, since 9 does not depend on the spatial variable: 



f{9{t)){9-9{t))d^l^ f{9{t)) / {9-9{t))di,^0. 
n Jn 



We also notice that for small 9 — 9 we have 



m^fi0) + n9)i9-9) + -f"i9)i9-9)' 
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where 9 lies between 9 and 9. We then have 
{II) 



f{9) [9 - 9{t)\ dti = / [f{9) - f{9{t))] {9 - 9{t))A^, 
n Jn 

= r(9) f [9-9{tfdf,+l f r{9)[9-m]' dfi. 

Jn ^ Jn 

Now we use the fact that 

f(9) = 1 - m(l + = 1 - m + m(l - m)9 + 0(f) 

which tends to (1 — m) as t — !■ oo uniformly in the space. Also 

f'(9) = m(l - m)(l + ^ m(l - to) 

uniformly in space as < — >■ oo. Putting these things together, we have obtained that 

-e{9{i)\ > TO / (1 + 9)-^-^ \^9\^ S-^™ da: - 2c[l - to + £:[6i(i)] 

> to[1 + e(t)r-^ [ I Vep 5*2™ dx - 2c[l - m + e(i)] f [6'(t)] 
Jn 

notice that in the limit to — > 1, the last term disappears, and we recover the standard L^— weighted 
estimates that hold in the linear case, see subsection 3.1 When m > 1 we obtain, using convergence in 
relative error as proved in |3]: 



(3.16) 



_d 
'dt' 



(3.17) 



dt 



£[9it)] > TO / (1 + 



ive^l^ S^"" dx - 2c[l - TO + eit)] £[9{t)] 



> to[1 - e{t)r-^ [ I Vep S-^" dx - 2c[l - TO + eit)] £[9{t)] 
Jn 



as claimed. The fact that the r.h.s. is positive for sufhciently large time follows again by convergence 
in relative error. □ 



3.3 Weighted Poincare Inequality and first rate of convergence 

In order to get a rate of decay for £[9{t)] we shall need a suitable version of the weighted Poincare 
inequality adapted to our problem, that we formulate next: 

GWPI: General Weighted Poincare inequality. Given to, il, there exists a constant K > 



such that for every stationary solution S > of problem (3.101 with constant c = 1, and for every 



9 e WQ'^{n,S^'''dx) we have 

(3.18) [ S^"'\\79\^dx>K [ S"'+^\9-9\^dx. 

Jn Jn 

Note that K depends only on to and Q. In order to apply this property to positive solutions Sc = Sm c 



of the elliptic problem (3.101 with constant c 7^ 1 we use the transformation 
(3.19) Sc{x)^ ^iSl{x) 

that produces a solution 5*1 of problem (3.10) if /i = c^^'^"™). Therefore, we have 



/ Sf\S/9Ydx^^i''^ / Sf\S/9Ydx>K^i^'^ / Si+"'\9~9Ydx^K^l'^-^ / S';^+^\9 - 9Y dx . 
'n Jn Jn Jn 
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In conclusion, the GWPI is formulated for problem (3.10) with constant c 7^ 1 as 

(3.20) / sl'''\ve\^dx>Kc f s'j'+^ie-efdx. 

Jn Jn 
Main assumption. We now make an assumption that is crucial for the rest of the paper 

(3.21) i^m - 2(1 - to) > Ao > 0. 

The rest of the paper will be based on deriving the consequences of this assumption and on the other 
hand, on justifying that under suitable conditions the assumption holds. As a first hint in the latter 
direction, if we consider the linear case and compare formula (3.8) in Appendix[6j3, with formula (3.20) 
we see that the latter holds with m = 1, c = Ai and K = (A2 — Ai)/Ai, hence Aq = if > in (3.21), 
and finally Aoc = A2 — Ai. When we deal with the PME case, this assumption is always satisfied since 
TO > 1, hence we do not care about the expression of the constant K in the GWPI. 

Here is a first important consequence of that assumption in the FDE case: 



Theorem 3.3 Assume that (3.21 ) holds for a given m G (to^, 1) and given ft C M.'^ . Let he a global 
hounded and positive solution to equation (3.11). Then for every 7 < 70 AqC there exists a time 
to > such that for all t > to 



(3.22) 



£[d(t)]<£[d(to)]e-''('-'''l 



Proof. It is immediate after the derivation so far. Assuming now that (3.211 holds we go back to 
formula (3.16) to get 



--£[eit)]>c{Km[l + e{t)Y--' 



2[i-TO + £(<)]) £[e{t)]. 



As e{t) — when t — )• 00, the conclusion holds. 
Next we deal with the simpler PME case where we just need a GWPI with any constant to get a rate. 



Theorem 3.4 Let to > 1 and let 9 he a global hounded and positive solution to equation (3.11 1. Then, 
for all (5 < 2 + ^"^^ there exists a time ti depending on to, d, j3 and on the constant > of the GWPI, 
such that 



(3.23) 



£[e{t)] < £[0{ti)] e-^(*-*i) for all t > t^. 



Proof. It is immediate after the derivation so far. Recall that c = 1/(to — 1) > in this case. We go 
back to formula (3.12) to get, for any 6 > Q and any t>tg large enough: 

d 



1 



At 



£\d{t)] > (Km\l - e(t)r-^ + 2[to - 1 - e(t]]) £\e(t)] 

m — 1 ' 



1 



- {Km\\ - e(t)]™-i - 2e(t) + 2[m - 1]] £\B{t)\ 



> 



Km 

TO — 1 



-(5 



£[0{t)]. □ 



3.4 Norm decay 

There is one more step to perform, since the entropy decay does not automatically imply the decay of 
the weighted L^-norm, because the mean value 9 is not constant along the nonlinear evolution under 
consideration hence it must be also controlled. This is also related to the fact that we can not specify 
boundary conditions for the relative error function 9. We first deal with the main case TOs < to < 1. 
Recall that S = Sc,m- 



18 



Theorem 3.5 Under the assumptions of Theorem 3.3 Then both entropy and the L^~norm, decay 
exponentially with the same rate 7 < 70 = AqC. More precisely there exists a constant k > such that 

(3.24) [ \e{t)\'^S^+"'dx < K£[e{t)] < Ke-'^'^'''^^£[e{ti)] 

Jn 

for all t > ^1 3> 1, where k depends on to, d and £[0{ti)], cf. the end of the proof. 

Proof. We first observe that since ||0(t)||oo ^ as t 00, hence also 9{t) ^ a,s t ^ 00, we can always 
assume \9\ < 1/2 and \9\ < 1/2, for all t > to. Next, we deduce the differential equation for 



m = 



^6l(t,a;)S'i+™dx 



X 



using the equation Ot = S'~('"+i)V • {S^"'V{1 + 9^) +cf{9): 

9'{t)^ , ^ r oii . 4: / e{t,x)S'''+Ux 

In ^ da; dr ^ /^^ 51+™ da; dt Jn 

= r . / V •(^""VCl + da; +^—^— / f{9)S"'+^dx 

J,,S^+"-dxJn ^5i+™da;A/^ 

^ c/^/(0)5"'+idx 

J,,S^+^dx 

where f{9) = {1 + 9) — {1 + 9)™ . By convexity of / we have that / lies above its tangent at the origin, 

f(9)>il~m)9, 

so that 

o'it) = ' Yffl+^CI'^'^ - - - ' 

where in the first step we have used Jensen's inequality since / is convex. An integration over [s, t] C 
[to, 00) gives 

9{t) > 6'(s)e'^(i-™)(*-") 

which implies that 9{s) < for all s > to, otherwise we get a contradiction with the fact that 9{t) — > 
as t —>■ +00. Next by Taylor expansion we get 

(3.25) f{9) = f{9) + f'{9){o -~e) + \no){e - W 

with 9 = a9 + {I - a)9 for some a G (0, 1). It is easy to check that -1/2 < 6^ < 1/2 since \9\ < 1/2 and 
1^1 < 1/2, so that 

mi^-m)(iy ™ < f'i9) = ""^^r""^ < to(1 - to)22-" 

so that 

<cf{9)+ccp^^2^'-m 
~ 5^1+'" da; ^ ' 



since we recall that 



f {9){9~ 9)3^+"' dx^O and £[61] ^ / [^(t) - 0(t)|'5i+" dx . 
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Finally, we recall that by Proposition |3.3[ for t > to we have that 

£[d{t)] < e-^(*-*«)f [0(io)] < e-^(*-*°) , 

so that 

'7^(1 - 22-mg-7(t-to) 



f(9) + fcoe--^^*-*") 



fid) 



Now define the function z{t) by the relation f{z) + fcoe-'^^*"*") = 0, i.e. z{t) = /"^(-fcoe-T^*"*")) and 
z{t) — ^ as t — >■ oo; it turns out that 

z'{t) > = /(z) + fcoe~^(*~*''). 

Hence z{t) > 9{t) for all t > whenever z(io) > (^{to)- But this fact is in contrast with the fact that 
9{t) — i' as t — > oo: on the line {t, z{t)) we have 6*' = 0, so we define the regions 



Z+ = {{t,y)\z{t)<y<0} 



and 



Z- ={{t,y)\y<z{t)} 



so that 9' > on Z+ and 6*' < on Z^. As a consequence, if 0(s) e , for some s > to, then 
^(t) e for all t > s, since in Z_ we have 0' < 0, therefore 9{t) can not go to zero as t — )■ oo, which 
is a contradiction. Finally we have proved that 9{t) € Z~^ for all t > to, which is what we need to 
conclude that 

> 9{t) > z(t) = /-i(-A:oe^^('-*°)) 



that implies 

since for \s\ < 1/2 we have 



\9{t)\ < fcie-T(*-*°) 



<max{|(ri)'(0)|, |(ri)"(S)|^}|.|<fci5. 
Now we conclude by observing that for all t > to we have 

l f \9{t)\^ 3^+"" dx < f \9{t)-9{t)fs^+"'dx+ f \9{t)\^ 8^+''' dx ^ 2£[9] + \9it)\^ f 3^+"" 



dx 



< 2e-'^'^'-'°'>£[9{to)] + /cie-^'^^*-*"' / S^+'"' dx < fcae-'^^*"*") 



which concludes the proof. □ 

Next we deal with the case m > 1. Recall that S is the unique stationary state. 

Theorem 3.6 Under the assumptions of Theorem 3.4 Then both entropy and the h^ — norm, decay 
exponentially with the same rate, more precisely 

(3.26) / [^(Ol'^i+^dx < 2 ff[0(ii)]e-('5-2)(*-*i) + |(?o|' / 51+" dx) e-2(*-*i) 

Jfi \ Jn J 



for all t > ^ 1, as in Theorem \3'. 
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Proof. The first part of the proof is identical to the one of the previous Theorem 3.5 so that we arrive 
to the differential equation for 9 

-, , d L0(t,x)S^+"'dx c L /(6')5""+i da; 

^ ' dt 51+™ da; S^+'"' dx - ^ - v ; 

where f{9) = (1 + 9) — {1 + 9)™-, c — l/(m — 1), and in the second step we have used the concavity of / 
together with Jensen inequality, and in the last step we have used that / lies below its tangent at the 
origin 

f{e)<-{m-l)9. 

An integration over [s,t] C [ii,oo) gives 

9{t) <9{s)e-^^-''^ and \e[t)\ < p(s) je-^*"'') 
where ti is as in Theorem |3.4| Moreover we have that 



-^f^o(f) o(f)\' ^'^"'^^ _ Ia\e{t)\'s^+"^dx 1 



|2 



/^^i+'-da; 2 7^1^^ ^^1/^,^1+™ da; 2j^S^+-^dx 2 1 
Combining the above result with the entropy decay of Theorem |3.4| 

£[9{t)] < £[9{ti)]e-^'^^-*'^ for all t > h. 

we obtain 

\9{t)?S^+'^dx^2E[9{t)] + 2\9{t)? [ 8^+"" dx 



n 



<2(^£[0(ii)]e-(^-2)(*-*i) + po|'y 5i+™dx^e-2(*-*i). g 



4 Stationary solutions and their limit as j9 — > 1. 

Let 1 < p < Ps and let Up be a solution to the elliptic problem 

-AU = Ap f/P in n 

(4.1) <( t/ > in f7 

u = on arj 

where Ap > if 1 < p < and Ap = Ai for p = 1. We are interested in the relation between solutions 
of the elliptic equation for different values of p € [l,ps), in particular we would like to see whether the 
limit V limp^i Up exists and under which conditions it is the ground state of the Dirichlet Laplacian 
$1 on 17. The existence of a limit depends on a normalization that we will discuss below. 

It is well understood by subcritical semilincar theory that positive weak solutions of the above elliptic 
problem are indeed classical solutions up to the boundary. Weak solutions can be defined as follows: a 
function Up G W^Q^'^(r2) is a weak solution to the elliptic problem (4.1 ) if and only if 



(4.2) [ [VC/p • V(p - XpUPif] dx = 

Jn 



for all (p e Wq'^{Q,). Notice that when p — 1 there is a positive solution, unique up to a multiplicative 
constant, while when p > I uniqueness is not always true, it depends on the geometry of the domain. 
The difhculty in understanding the limit oi Up as p — >■ 1"*", relies indeed in the lack of uniqueness and 
on a scaling property typical of the nonlinear problem. In the case of uniqueness, for example in the 
case when f2 is a ball, solutions are variational, in the sense that they are minima of a the functional 
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||V/7||2 under the restriction ||?/||p+i = 1, but when the uniqueness is not guaranteed, solutions are just 
critical points of such functional. 

One can also easily see that the constant Ap > in the nonlinear problem can be manipulated by 
rescaling, because if J7p is a solution with parameter Xp^(i), then C/p^(2)(a^) = ^^Z^^""'^-' J7p (2)(x) is a 

solution with parameter Ap.(2) — /iAp.(i). In any normed space ||?7p,(2)|| = A^^^*-^ "'^''||C^p,(i)||- This means 
that scaling allows to fix the norm of a solution: changing the norm by a factor /i^/'?'^^) by scaling is 
equivalent to changing Ap in the equation by a factor /i^^. 

Assumption throughout this section. Let us fix Ap as the factor for which ||?7p||p+i = 1, so that, 
using Up as test function, we obtain the following identity 

(4.3) ||VC/p||2 = Ap||;7p||^t} = Ap, 

so that it is equivalent to prove that Ap — > Ai or to prove that ||Vt/p||2 — > ||V<i>i||2, when p — > 1. Recall 
that $1 has unit L^-norm. 

We state now the main result of this section. 



Theorem 4.1 Let Up be a family of solutions of Problem 4.1 with p G [l,ps), \\Up\\p^i = 1 and let 
Xp > be chosen according to (4.3). Then as p ^ 1, Ap — >■ Ai, Up $i in L°°(r2) , Vt/p —> V$i in 

(L^(ri))''. Besides, there exist two explicit constants < cq < Ci such that 

(4.4) cg-^Ai < Ap<c?-^Ai. 

Moreover, there exists constants < fco(p) ^ ^i(p) such that ki{p) ^ 1 as p — > such that 

(4.5) ko{p) < < ki{p), for all x . 



4.1 Proof of Theorem 14.11 

The proof of the above Theorem will be divided into several steps. Our first result in this connection 
is the following. 



Lemma 4.2 Let Up be a solution of Problem 4.1 with p G [l,ps), ||J7p||p+i = 1 and let Xp > be chosen 
according to (4.3). // there is a constant A> such that 

< Xp < A < oo , 

then Up $1 in L'(f7) for any < q < 2* , and Ap — !■ Ai. 



Proof. Since Up is a solution to the elliptic Problem 4.1 with HL'pllp+i = 1 we have that the hypotheses 
together with the energy identity (4.3) give 

< Ap = ||V[/p||^ <A<oo 

which proves that Ap = ||VC/p||2 is uniformly bounded for all p G [l,ps). Hence we can guarantee 
that there exists a subsequence VUp^ that converges weakly in L^(r2) to a function W . Moreover, by 
Kondrachov's compactness theorem there is a (maybe different) subsequence Up^, that converges to V 
strongly in any L'^ with 1 < g < 2*. Strong convergence implies that |1F||2 — 1, hence V can not be 
identically zero. Moreover, it is well known that in this case W = W. Next, we show that U^^ — >■ V 
inii(17): 



|C/P; -V\dx< 



H:-Up^ 



dx 



\U, 



Pn. 



V\ dx 



1 1 Up^ dx - 



\U, 



Pn 



V\ dx 
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where the second integral converges to zero since Up^ — > V in L^(r2), while for the first we have to use 
a numerical inequality: 



II < 



a" 1 
b a 



for all 



a>0, 0< z <b 



which we will prove at the end of the proof. Using the above numerical inequality in the first step for 
a = J7p„ > and < z = p„ - 1 < 6 = 2* - 2 = 4/(d - 2), we obtain 



|C/;;-i -l|C/p„ da;<(p„-l 

1 



b 



1 



< 



< 



Pn 

b 



b\n\ + 



Up^ dx < 



Pn 

b 

b+l -1 



1 



u: 



b+l 



b) dx 



2*\n\ + \n\ 



\VUpJ'^dx 



Pn 



1 r 



2*|r!| + |r!|i-^ (52Ap)^ 



< 



2*\n\ + \n\'-"-i^ {S2A)'^ 



since b + l<2* — l=ps — {d + 2)/(d — 2), we can use Sobolev and Holder inequalities, and the fact 
that the L^-norm of the gradient, or equivalently Ap, is uniformly bounded by A. 

Now, we identify the limits. For all 1 < p < pg-. 

(4.6) \p f UP^idx = - j ^i/^Updx^ [ Vf/p • V$i dx = - / C/pA$i dx = Ai / Up^idx 
Jn Jn Jn Jn Jn 

such equalities hold by the weak form of the equation satisfied by Up and since both Up and $i are in 

Taking limits 



, A„ 



Wq' (Q). Take any subsequence of p„ such that Ap^ converges. Let A = limfe^oo '-pn^ 
as fc — >■ cx) so that — > 1 in the above expression to get 

A / V^idx = Xi I V<^idx 



since we know that both Up^^ and Up^^ converge to V. We conclude that limfe_j.oo Ap^^ = Ai. Since 
this holds for all subsequences, this means that lim„_>.oo Ap„ — Ai. We are now ready to identify 
V = lim„_j.oo C/p„ , indeed we just notice that 



WV -Wipdx- Xi / Vipdx 



lim 



/ VL/p„ • V-^dx - Ap„ / UP;ydx 
Jn Ja 







where all the quantities have been shown to converge in such a way. The latter equality identifies V as 
the unique ground state $1 such that -A<I>i = Ai$i and ||F||2 = W^ih = 1- 

Now we prove that for any subsequence p„ — ?■ 1, we have Ap„ — > Ai, Up^ — !■ $1 and consequently 
Up" ^ $1. Suppose that there exists a sequence p„ such that lim„ Ap„ — A ^ Xi. We can repeat the 
first steps to conclude that there is a subsequence Up^^ that converges strongly in to some V and 
such that VUp^^ converges weakly in to 14^ = W . Moreover, we also have that U^^ ~^ ^ ■ Using 
formula (4.6), we have that 

V, / Ulll $1 dx = Ai I Up^^ $1 dx 
'' Jn ^ Jn 



and taking the limit as /c — cx) we get a contradiction, since limfe Ap 
once we prove the numerical inequality 



1 < max 



a" - 1 1 
b ' a 



< 



1 



: Ai. The proof is concluded 
for aU a > , < z < b. 
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we prove it first for a > 1: since the function f{z) — 1 is convex, it lies below the secant (a^ — 1)-^/^ 
for all < z < 5, hence the inequality ~ 1| — — 1 < (a^ — l)z/b when a > 1 and < z < b. 
When < a < 1, we see that — 1| = 1 — a^. We see that f(z) ~ 1 — = z/a when z = 0, and that 
f'{z) = — log(a)a^ < log(l/a) < 1/a, hence the desired inequality is valid also when < a < 1. □ 

In view of this result, we need to prove the upper bound Ap < A. The proof is based on an idea of 
Brezis and Turner [TS]. This method relies on an Hardy-type inequality, which holds for a large class 
of domains, but only in the range 1 < p < {d + l)/{d — 1). If one wants to deal with the full range of 
exponents 1 < p < Ps, one has to proceed as Gidas-Ni-Nirenberg [21] when the domain is convex, or as 
DeFigueredo-Lions-Nussbaum [20 which extend the ideas of |28] to more general domains. 

Proposition 4.3 The following Hardy-type inequality holds true whenever Q has a finite inradius and 
satisfies a uniform exterior ball condition 



(4.7) 



< Hr,d\\^f\\ 



iffew^'^in), o<q< 



2d 



d-2 + 2r 



and < r < 1 



where $i is the unique positive ground state of the Dirichlet Laplacian on il, and H^^d is a suitable 
positive constant that depends only on r,d and \ and is given at the end of the proof. 

Proof. The proof is obtained by combining the standard Hardy inequality 



/ 

dist(-,af7) 



<i7o||V/||2 ioT all f eW^'\n), 



which holds whenever has a finite inradius and satisfies a uniform exterior ball condition, see for 
instance Section 1.5 of [1^. We combine such Hardy inequality with the standard Sobolev imbedding 
II/II2. < 5|||V/||2 as follows: 



dist(a;, dflYq 



dx 



f 



rq 



< 



da; 



--{a) 
'-(b) 



^ Vdist(a;,917)''« 
dist(.T,9rj)2 



da; 



<(.) [How^fhr \n\'-'''^sr^>''\\\7f\\t^^'' 



dx 


7 


[/ 


rq 
2 


[J 


f 2 

T 


-2+2r 


,(1- 
'2 


2d 





(l-r)9 



da; 



where in (a) we have set 7 — 2/rq > 1 and 7/(7 — 1) = 2/(2 — rq), for any < r < 1, while in (b) 
we have used the above Hardy inequality and we have estimated the second integral with the Holder 
inequality 



2(l-r)<j 

/ da; 



= 11/11^^ <N^^-'^^'(^--)||/lir^^' 

2-rq 



and in order to use Sobolev inequality, we need 
(1 — r)<?7 2(1 — r)q 



7-1 



< 2* 



that 



< 



2d 



rq 



We have obtained so far 
(4.8) 



dist(a;, dfl) 



rq 



■ dx 



1 d-2 + 2r 



d - 2 + 2r 
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We conclude the proof by noticing that, since there exists two positive constants cg, ci depending only 
on the dimension d > 3, such that 

codist(-,af7) < $i(-) < cidist(-,5f}) in fl. 



We combine the above lower bound together with (4.81 to get the desired inequality (4.71 

1 1 



$1 



rq 



dx 



< 



1 



dist(x, di}) 



rq 



dx 



H, 



We are now ready to prove the upper bounds for Xp. 



Proposition 4.4 Let I < p < {d+ l)/{d— 1) and d>3 and \p be such that \\Up\\p+i ~ I, as in (4.3 1. 
Then the following upper bound holds true 



(4.9) 



[d + l-p(d-l)l(p+l) _2p_ 



Proof. Testing the equation (4.1 1 with $i yields, as in (4.6) 



Xp [ C/;$i dx = Ai / Up^i dx or ^ [ f/;$i dx^ [ [/p$i dx 
Jn Jn Ai Jq 



that gives 



Ai 



C/P$i dx = / Up<t>i dx< [ C/P$i da; 



$1 dx 



where we have used Holder inequality in the last step. We have obtained that 



(4.10) 

Next we calculate 



UP<^>i dx < 



Xr. 



$1 da; . 



1^ / UP+'dx= f {U^^^) 
Jn Jn 



da; < 



(a) 



f/P$i da; 



n ^1 



< 



'-{d) 



Xr 



$1 da; 



$1 dx 



Up 



g(l-Q) 



< 



AiA- 
A,] 



<i>i da; 



(i/,,rf||v;7p||2) 



g(l-a) 



(yHr,d^p ) 



(1-q)p+1 



where in (a) we have used Holder inequality with conjugate exponents 1/a and a/(l — a), for some 
a G (0, 1) to be fixed later. We have also put 

p + 1 = ap + (1 — a)q and a = r(l — a)q 

which is equivalent to fix the values of q and r (as functions oil < p < Ps and < a < 1) as follows 

1 , 1 a 



q^p + 



I- a 



and r 



{l-a)q l + {l~a)p 



we notice that < r < 1 since < a < 1 and p > 1. In (6) we have used the upper bounds (4.10 ) , while 
in (c) we have used the Hardy inequality (4.7), for which we have to check that q < 2d/{d — 2 + 2r), 
that is equivalent to 



< 



2d 



I- a ~ d-2 



l + (l-a)3: 



d+2-{d-2)p 

that is a< — , , , ^ < 1 

- 2{d+l) - {d-2)p 
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since p < ps = {d + - 2) and hence < d-\- 2 — {d — 2)p < 2{d + 1) - (c? - 2)p. Finally, in the last 
step {d) we have used the identity ||V[/p||2 = ^p||^p||p+i = and the fact that q{l~a) = (1 — ct)p+l. 

We have obtained the desired upper bound (4.91 



pg _ (1 — q)p + 1 ay 



for any a € (0, 1) such that 



pa 



(1 -a)p+l 



> 



that is 



a > 



p-1 



that is 



p-1 d + 2-(d-2)p 
< a < 



p 2{d + 1) - {d - 2)p 

which is non empty only when p < {d + — !)• Letting now a = 2/{d + 1) gives the desired upper 
bound (4.9), since 



pa 
p^ 



{l-a)p+l _ -{d - l)p2 + 2p + d+l _ [d+l-p{d-l)]{p+l) 



2(d+l)(p-l) 



2(rf+l)(p-l) 



□ 



At this point we are able to prove the first part of Theorem 4.1 



Proposition 4.5 Let Up be a family of solution of Problem 4.1 with p G [l,Ps), — 1 and let 

Xp > be chosen according to (4.3). Then as p 1, Ap — > Ai, Up — ?> $i in L°°(r2) , Vt/p — > V$i in 
(L^(f2))'^. Moreover, there exist two explicit constants < cq < ci such that 



(4.11) 



^Ai < Ap < ^Ai , 



Proof. The only thing that remains to prove is the convergence in L°° and the convergence of the 
gradients. We know that Up and $i are Holder continuous on the whole fi, and the C"-norm of both 
functions is uniformly bounded, say ||t/p||cc.(f2-) + ||$i||c°(f2) ^ Moreover combining the upper 
estimates for Ap of Proposition 4.4 (or of Proposition 4.7 when there is uniqueness of the stationary 
state) with Lemma 4.2 gives the convergence of Up — >■ <i>i in L*(r2) for any 1 < q < 2*. By the 

we get that Up — $1 in L°°(f2), since when p — )■ 1 we have 



interpolation Lemma 



5.9 



|f/p-$l||co <C\\Up-^, 



l-i5 I 
C°(0)l 



Up-'i'iWt <CK'-'\\Up^^if^^O. 



As for the gradients, we use the Hardy inequality (4.7) with r — 1 and q — 2 applied to J7p — $1 G 
W^^\n), to get 



dx<ffi,2 / \V{Up~^i)\ dx 



\Up-^i\ 
$2 

Next, we analyze the right-hand side: 

[ |V(;7p - $i)f dx = [ \\/Up\^dx+ [ |V$ifdx-2 [ Vf7p • V$i dx 
Jn Jn Jn Jn 

(a) - Ap /" [/P+i dx + Ai ( $2 da; - Ap / U^ $1 dx - Ai / Up $1 da; 
Jo, Jo, Jn Jn 

UP\Up~^i\dx + Xi [ |C/p-$i|$ida; 



< A, 



(6) < Ap 



[/P+i dx 



_, P 

P + 1 



da; 



+ Ai 



da; 



$t da; 
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where in (a) we have used formula (4.6) in the form 



2 [ \/Up ■ V$i dx^Xp [ U'P $1 da; + Ai ( Up $i dx . 

As for (6) we have used Holder inequality with conjugate exponents ij) + l)/p and p + 1 for the first 
term, and Cauchy-Schwartz inequality for the second term. Putting all the pieces together we have 
obtained 



$1 



< K\Wp\\l+i\\Up ~ <i>i||p+i + Ai||c/p - $i||2||$i| 



p-i. 1+ 







since we already know by Lemma 4.2 that Up — )■ $i in L'(il) for any 1 < g < 2*, and we also know 
that p + 1 < 2* since p < Ps- The sole requirement of Lemma 4.2 is that Xp < A, and the uniform 
upper bound for Xp is guaranteed by Proposition 4.4 when 1 < p < {d + l)/{d — 1) or by Proposition 



4.7 when Up is a variational solution and 1 < p < Ps- □ 



The last step in proving Theorem |4. 1 1 consists in comparing solutions corresponding to different p and 
Xp, more precisely to show that there exist constants < ko{p) < ki{p) such that ki{p) — 1 as p — >■ 1+. 
$1 is the corresponding ground state, towards to Up converges as p — > 1. 

Proposition 4.6 Under the running assumptions on Up and $i, there exist constants < fco(p) < 
ki{p) such that ki{p) 1 as p 1+, such that 



(4.12) 



r , ^ U^ix) 



<ki{p), 



for all X fl . 



Proof. The proof is divided in several steps. 

• Step 1. Convergence of the quotient in an inner region. Proposition 4.5 implies that Up/^i — )• 1 in 
any inner region in which $i > ct > 0. In the sequel we will construct a special region as follows. By 
Lemma [2.41 we know that there exists a ^ > such that 



Vs = {xen : d{x) < 6} 

such that d{x) = dist(x,9J7) S C^(r2^) and d(x) is Lipschitz with constant 1, i.e. \d{x) — d{y)\ < \x — y\, 
and < c < |Vc?(a;)| < 1, ~K < Ad{x) < K in ^Ig. In the complement, ilg = ft — Vs, we know that 

o- < Up{x),<^i{x) < M for all x e ng, 

so that, given e > for p sufficiently close to 1 we have 

(1 - e)$i(a;) < Up{x) < (1 +e)$i(a;) for ah x e fig . 

It remains to prove that the above inequality extends to the thin region Vg — ~ fig, and this will be 
done in the next steps. 



• Step 2. Upper comparison near the boundary. The upper estimate for Xp of Theorem 4.1 reads 
Xp ^ Cj' Ai. Since we are working in the thin annular domain Vg = D,\Qg, and we know that $i = 
on dfl, we can assume (eventually by taking a smaller (5 > 0) that (ci(l + £)$i)^^^ < 1 in Vg. As a 
consequence, we have that 

Ap((l + < Ai(ci(l + e)<i>iY'\l + < Ai(l + e)$i . 

This allows to compare Up and $2 = (1 + e)$i on the thin set Vg. We have < Up ,{1 + e)$i < M in 
Vg. The respective equations are 



-AUp = Xp m , 



in Vg , 



-A$2 = -A(l + e)$i = Ai(l + e)$i > Ap((l + £)$i)^ = Ap$^ in Vg , 
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and the boundary data 

$2 C/p = , ondn, 
$2 > C/p > , on dfls . 



We want to apply Theorem 6.5 to obtain the comparison Up < (l + e)^i — $2 in Vg. We need a the 



following smallness condition on Vs: 

The above condition can be fulfilled just by choosing S sufficiently small, and we can always do that 
independently of e small. 

• Step 3. Lower comparison near the boundary. This part consists of two comparison arguments. First 
we observe that we can compare Up with a suitable harmonic function on Vg-, namely 

-AC/p = Ap f/P > , inVs, 
-AC/ = in Vs , 

and the boundary data 

u = Up^o, on an, 



[/ = (1 - e)$i <Up, on OQs ■ 
we apply standard comparison to get U < C/p in V5. 
Next we want to prove that U > {I — 2e)<i>i on Vg if S is small enough. 
We define the function w — {I — £)<&! — U that satisfies 

-Aw = Ai(l -£)$!< co(5 , inVs, 
w — , on dVs , 

and we compare li = w/^cqS) with W, which solves the following problem on the whole = il^ U V^: 

-AW = 1 , in , 
W ^0, ondn. 

By comparison, we have that w < W in Vs, which means 'w{x) < coSW{x) . Moreover, we know that 
the function W satisfies W{x) < c\ d{x) < C2^i{x) since we know that $1 > cd{x). Summing up we 
have proved that 

w{x) < coSW{x) < ci6d{x) < C2S^i{x) , 
recalling that w ~ {1 — e)^i — U, the above inequality gives 

[(1-e) -C25]$i < U 

Now putting C2S < e we get the result, when 6 is small enough. 

• Conclusion. The two above steps imply that given e > there exist a (5 > and > I such that 
the above two steps hold and 

(1 - 2e)$i < C/p < (1 + £)$i in f7 . □ 



Combining Proposition 4.5 with Proposition 4.6 we get the full statement of Theorem 4.1 [] 



4.2 Additional bounds on Ap 

We shall also prove suitable lower bounds for Ap, both for the sake of completeness and because they 
will be used in Section [5.3[ These bounds are easier to obtain than the upper bounds. 
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(i) Using Up as test function, we obtain the global energy equality Ap||[/p||^^J = l|V[/p||2, that combined 
with the Sobolev inequality 

II/I!p+i < m^^-^Wfh' < \n\^^-^S2\\vf\\2 

gives, recalling that we have chosen Ap in such a way that ||C/p||p+i — 1, 



i \\Up\\l+i ^ 



|r2|p+i 2* \n\p+i 2^ 

We can rewrite the lower bound as follows 



U; dx 

n 



<5||lVC/p|l2 = 5|Ap|lC/p|lP+}=5|Ap 



(4.13) \ ^ < A„ andforn^'l ^ < Ai . 

(ii) Other lower bounds can be obtained by combining Holder, Poincare and Sobolev inequalities: 

liC/pg+r < \K\\lf\\Up\\r-'^ < (^||VC/p||^ w.th ^ . 1^ 

which gives 

(4.14) Ap - / |Vf/pp dx > -^\\Up\\l+, = Ai {X.Siy^ 

Jn (Aio|) 

since we have chosen Ap in such a way that ||[/p|jp+i — 1. 

The case of variational solutions. Other estimates for Ap can be easily obtained in the case in 
which solutions are minima of a suitable functional, this happens for instance in the case of domains 
fl for which the solution is unique, hence they are minima, since a solution which is a minima always 
exists as a consequence of Kondrachov's compactness theorem. 

When the solution of the Elliptic problem |4.1| are minima of a suitable functional, namely when we 
consider the homogeneous functional 



defined on Wg^'^(ri), and we seek for its minimum under the restriction ||7i||p+i — 1, we can define 
Ap = inf J.p[u] = inf / \Vu\^dx where Xp = \ue W^'^in) I ||ii||p+i = l) . 



Let Up € Xp be a solution to the elliptic problem 4.1 with Ap defined as above. Estimates in this case 
are simpler and hold for any 1 < p < Ps- 

Proposition 4.7 Under the above assumptions, if Up is a minimum for the functional Jp on the set Xp, 
then it is a positive weak (hence classical) solution to the elliptic Problem 4.1, Moreover the following 
estimates hold 

inf Jp [u] 

(4.15) ^s2X,r^<^ = --^^^<\n\'^ 

Ai mt Ji\u\ 
ueXi 

where Ai is the first eigenvalue of the Dirichlet Laplacian on Q, and S2 is the constant on the Sobolev 



imbedding from WQ^'^(r2). As a consequence, Ap — > Ai as p — > 1 



+ 
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Proof. It is a standard fact in calculus of variations to see that a minimum of Jp is a weak solution to 
the elliptic problem under consideration. We can now prove the upper estimate: 

Ap= inf / \\/u\^dx= inf ^ < ^^^^-^ ^<Ai|f7|FFT 

if we moreover assume ||^i||2 — 1 (not restrictive). We have just used the fact that A$i — Ai$i 

p— 1 

together with Holder inequality ||$i||2 < l^l^ll'J'illp+i- The lower estimates are exactly the same as 
(4.141 and we do not repeat the proof here. 



5 Convergence with rates for m near one 

The idea of this section is simple, but the technical details lengthy. Since the lower bound for the decay 
rates in the linear case m = 1 is just Aqc = A2 — Ai > 0, it must be also positive for m near 1 by 
continuity. Putting the details into this program is not so easy and we present below the part that we 
have been able to prove. The section starts by proving a suitable Poincare inequality. It continues by 
estimating the constant c that enters the elliptic problem and proving that it tends to Ai as it should. 
These two ingredients allow to state and prove our main results about convergence with rate. 



5.1 Weighted Poincare Inequality 

We recall that putting p = 1 /m and Sm = we get a solution Sm to 

Sm > in 
5'rr, = on dO. 



for which we know, by Theorem 4.1 that 1 = = ||<I'i||2 and, as TTi y 1 , \m — ^ Ai, and 



kQ{l/m)<^i < 5™ < fco(l/m)<I>i, with lim„j^i ki{l/m) = 1. Hence setting ki{m) := fc,;(l/m) we have: 

(Hm) For any rus = {d — 2)/(d + 2) < m < 1 there exists constants, ki{m) with lim„j-|-i ki{m) = 1, such 
that the stationary solutions Sm{x) satisfy the bound 

(5.1) fco(m)$i(2:) < S';'J(x) < fci(m)$i(x) for any a; e H . 

Theorem 5.1 (Weighted Poincare Inequahty) Let f G WQ''^{il) and g = Let Sm be a 

weight satisfying (Hm). Then the following inequality holds 

(5.2) , , ^tii r 111--^ < / \Vgf Si- dx 
fci(TO)-^||S'„||oo Jn Jn 



where A = A2 — Ai > is the optimal constant in the intrinsic Poincare inequality (3.8 I ana 

j^gSl^-dx 



9 



In " '^^ 



Proof. Notice first that, by (Hm), ^l{x) < S"^ /ko{m)'^. As a consequence, 

(5.3) / \Vg\^ dx < —1^ / I dx . 

Jn ko{my Jq 
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Moreover 
(5.4) 



o2m 



l+m 



1 1— m 

oo 



fcl(m)2 - fci(TO)2||5„.|| 

where we used again (H^) and the fact that S^{x) > ||'S'm||^~^5'^^™ , vahd since m < 1. Therefore 



\9 - g^A' S'J"' dx 



(5.5) 



/ |g - o,!,, 1^ $?da; > / 

ki[my\\Sm\\oo Jn 



where the last inequahty follows by Lemma [5?2| 



Putting together the latter inequalities yields (5.2) 

We now recall the following well-known result which has been used in the above proof. Hereafter, 
//J := j-^ f dfj, where /i is any nonnegative bounded measure. 

Lemma 5.2 Let f G L^(X, d/i), with fi{X) < oo. Then we have 

(5-6) 11/ - Uk^xA,) ^ 11/ - ^llL=(x,d^) : for allceR. 

Proof. By expanding the square 



11/ - c|lL^(x,d,,) = J^\f- ^I'd^ = ^(^) ^ - 2c/ + c2] 



= t^{X) 



f 

X KX) 



/llL2(X.dAi) 

Kx) 



WfWvHxAt.} - t^ix){f^y = 11/ - fXmx.df.) ■ □ 



The version we will really need is a variation in the use of different weights. Let S'c,m be a positive 
solution of — AS*™ — cS on Q, vanishing at the boundary, obtained as asymptotic profile for an evolution 
with fixed initial data uq and variable m < 1. 

Theorem 5.3 Let / e Wg^'^(ri) and g = //*&!. Let Sc.m be as above. Then the following inequality 
holds 



(5.7) 



fci(m)2||S,,, 



I \g-g?Slj:^dx< f \Wg\'s',Z<ix 
loo Jo Jn 



where A — X2 ~ Xi > is the optimal constant in the intrinsic Poincare inequality (3.8 1 and 

In ff'^c.m™ 



g = 



In •S'c.m™ dx 



and we know that limm-^i ki{m) /ko^m) — 1. 
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Proof. We use the preceding result and the following observation: it easy to check the relations between 
Sin and Sc ni' 



Q — 



An 



and 



"f II 1— m 
~c,m\\q 



An 



Sn 



We now plug the above equalities in the weighted Poincare (5.2 1 to get 



Kko{my 



ki{mY^\\S,,m\\ 



1 — m 



Iff -31 



iv.9r 



that is exactly (5.7) since the factors c/A™ simphfy. □ 



5.2 Estimating the extinction time and the constant c 

We need to estimate the extinction time T — T{m, d, uq) from above and from below, to obtain bounds 
on the constant c = 1/(1 — m)T appearing in the rescaled equation ( 1.3 1 and the elliptic equation ( 1.5 ). 

Proposition 5.4 Let < m < 1 and u be the solution to the original Problem |1.1| corresponding to an 
initial datum Uq £ {Q) with r > 1 and r > rc = d{\ — m) /2. Then its extinction time T = T{m,d,Uo) 
satisfies the bounds 



1 [J^^Uoix)^,ix)dxY ^ ^ {r + m-lf (Ai5|)' 

Ai [/^$,(a;)da;]'~" " 4m(r - 1) Ai 

Taking r = 1 + m, which amounts to ask m > in^ = {d + 2)/{d — 2), we get 



(5. 



(5.9) 



1 [J^^uo{x)<i>i{x)dxY 

1 1 — m 



< (1 - m)T < 



(Ai5|) 
Ai 



Ai [^a>i(x)dx] 
Corollary 5.5 If uq E L'™+^(ri) and m > rUg we have 
(5.10) lim {l-m)T{m,d,uo) 

m— f 1 

hence c — > Ai as m — > 1. Moreover, (c/XiY/'^^^'^^ = 0{l) as to — > 1 



^*olll+m ■ 



1 

A^ 



Proof of the Proposition. We begin with the lower bound. We take $i as test function and consider the 
solution u to the original Problem |1.1| 

u^ = A(w") in(0,+oo)xl7 
u(0, x) — uo(x) in i7 
u{t, a;) = for T > and x € dil 



and derive the integral 
d 



dr 



u(t, a;)$i(x) da; 



(Au'"(T,a;))$i(x)da; = 
Ai / u"'{T,x)<^i{x)dx < Ai 



/ M'"(T,a;)A$i(x)da; 
Jn 

u{t, x)^i{x) dx 



m"(t, a;)Ai$i(a;) da; 

1 — m 



$i(a;) da; 



where we have integrated by parts since both u and 3>i are zero at d^l, and we recall that Ai$i = — A<i>i. 
Integrating the differential inequality gives 



u{t, a;)<i>i(a;) da; 



- 1— m 




< 


[/ 







u{s, a;)$i(a;) da: 



+ Ai(l - m) 



-| 1 — m 



$i(x) dx 



\t-s\ 
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for any < s,t < T. Letting s — T and t = gives the lower bound (5.8). 



The upper bound follows by using Holder, Sobolev and Poincare inequalities in the form 



11/11? <ll/ll^?ll/llf-''^< 



^'^''^^^ !|v/!|^ with 



Ai 



2q 



and 2<q<2* 



We are going to use the above inequality for q = (2r)/(r + m — 1) and f = u 2 ^ noticing that 
q e [2, 2*] if and only ii r > Vc ~ d{l — m)/2, which is 



Ai 



d(l-n.) 

(Ai5|)^ 



Mr 



r+7n— 1 



< 



Differentiation of the L*"— norm gives when r > 1 and r > = — m)/2 



(5.11) 
d 
dt 



4mr(7' — 1) 
(r + m — 1)2 



(r + m — 1) 



(Ai5|)^ 



which is a closed differential inequality of the form Y' < ~KY^ ^, which integrated on [s,t] gives 



Hs)\\ 



Amr(r — 1) Ai 1 — m , 
^ 4 ... ^, , ft-g) 



(^ + ™-l)'(Ai5|) 
which gives the upper bound, just by letting s — Q and t = T. □ 



5.3 Statement of the main convergence result 

The next step consists in showing that these inequalities allow us to apply the decay results of Subsection 



3.3 For that we need to estimate in a clear way the constant before the left integral in (5.71 and obtain 
a lower bound with a constant K independent of the particular solution S. This is possible thanks to 
the following estimate 



\Sr. 



X 



A 



-< 1 1 1 — m 



^c,m II 00 

which is true since \\S, 



Q II 1 — m 
„ n'-'mlloo — „ , , i-i 
C C ^T+T, 



A. 



(52Ai)5 \n\ 



p = 1/to, reads 



miii+m = 1 t)y construction and since the lower bound (4.14), rewritten for 

A,: 



Ai 



> (52A1) 2(1 + '") 



This means that the Weighted Poincare Inequality mentioned in Subsection |3.3| holds in the form 

(5.12) Kc [ ST+'\9-9\'dx< f ^f"|Vgpdx, 

Jn Jn 



with 



K = K{jn) = 



(A2- Ai)fco(TO)' 



(52A1)' \n\ 



At this moment, we see that the necessary condition to obtain decay is then 

F(m) := mK{m) - 2(1 - m) > 0, 

in other words, ni > 2/(2 + K{m)). But since F{1) — (A2 — Ai)/Ai > 0, it follows that there is an 
TOj < 1 such that for all < m < 1 we have that F{m) > 0. Note that toj changes with the geometry 



33 



of the domain. It may be objected that toj is given in a very imphcit way. However, given an estimate 
of the form ki{m) / ko{m) < C when iric < m < 1, then 1 — mj can be expHcitly estimated from below 
in terms of C, Ai, A2, ^2 and \ We shah provide suitable explicit bounds for ki{m) in a forthcoming 



paper [TJ], see also Theorem 6.4 



We can now state the rescaled version of the asymptotic convergence result. The rate will involve the 
expression 



(5.13) 



7o(to) 



1 



(1 - m)T 



A2 
Ai 



1 



fco(TO 

fci(m) 



l2 r 



(52Ai)^ \n\ 



2(1 -m) 



> 



where mj < m < 1, the constants ki{m) — >■ 1 as m — >■ 1, cf. Theorem 4.1 (rewritten there with 
p = 1/m). We recall that by Proposition 5.4 the quantity (1 — m)T appearing in (5.13) can be 



explicitly estimated from above and below in terms of uq and moreover, by Corollary [5^5] we have that 
lim l/(l-TO)r = Ai. 



m— > 1 



Theorem 5.6 (Rates of convergence, rescaled version) Let max{mu,TOc} < m < 1. Let v he the 

rescaled solution corresponding to an initial datum uq as in Theorem \2.\\ which converge to its unique 
stationary profile S. Let 7 < 70 then for all t > with to sufficiently large, we have the following 
entropy decay formula: 



(5.14) 



£[0{t)] < e^T(*-*°)£:[0(io)] 



In other words, the weighted if^-norm decays with rate 7, more precisely there exists constants Ki > 
and a time to > such that 



(5.15) 



\v{t,x)- S{x)\^ S{x)"'-^dx 



v{t, x) 



Six) 



1 



S{xy+"'dx < Koe" 



-7(t-*o) 



Moreover for all q £ (0, 00] 

(5.16) \\yit,-)-Si-)\\g < «:ie-^(*-*«) 

where the constant ki depends on m,d and uq and on the uniform bounds on the C^ — norm ofu{tQ). 
Remark. The constant 70 satisfies 



lim 70 (to) = A2 - Ai , 

m— ^1+ 



as a consequence of Corollary |5.5| 



Theorem 5.7 (Rates of convergence, original variables) Let max{mu,TOc} < to < 1. Let u he 

the solution to Prohlem \Ll\ let T = T(to, d, uq) be its extinction time and letUx be as in Theorem 2.1 
so that u{t) /IAt{t) — > 1 as t ^ T . Then, for any 



7 < 



1 



A2 
Ai 



(1 — to) 

there exists a constant k > such that 

u{t, ■) 



- 1 



(5.17) 

or equivalently 
(5.18) 



U{t,-) 



koijnf 
ki(my 



(52Ai)3 |r!| 



- 2(1 - to) 



\u{t, x) - U{t, x)\^ S"'-^ dx < Ko 



T~T 

T 



T-T 



T 



r+7 
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for all to < T < T , where kq depends on m, d and Uq. Moreover we have that for all q G (0, oo] 

(5.19) \\u{T,x)-UiT,x)\\,<Ki 
where the constant ki depends on m,d and Uq and on the uniform bounds on the C^ — norm of u{to) 



T-T 
T 



We comment that the weighted convergence of (5.181 is somehow stronger than the non-weighted 
LP— norm convergence, since the weight S"^~^ is singular at the boundary. 

The main resuh for the Porous Medium Equation reads: 

Theorem 5.8 Let m > 1, let v be a the rescaled solution as in Subsection \2.4\ that converges to its 
unique stationary state S, and let 9 = v/S. Then, for all < f3 < 2 + ^zri there exists a time ti 
depending on m,d,l3 and on the constant K > of the GWPI, such that the entropy decays as 

(5.20) £[0{t)] < £[0{ti)] c^^(*-*i) for all t > ti. 
Moreover for all q G (0, oo] 

\\v{t,-)-S{-)h,^n) < «ie-(*-*«) 

for all t > ti 1, where the constant K2 depends on m, d and uq and on the uniform bounds on the 
C^ — norm of u{t()). In original variables we obtain that for all q G (0, oo] 

||u(r,-) -Z^(T,-)|lL<j(n) < ''ii 1 • 

(1+7-)!+™-! 

In order to conclude the proof of the above theorems, we need an interpolation Lemma due to Gagliardo 
gH, cf. also Nirenberg, (33 p. 126]. 

Lemma 5.9 Let \, fj, and v be such that —oo < X < fi < h' < oo. Then there exists a positive constant 
C\,^i,iy independent of f such that 

(5.21) 11/11'^/-' < Cx^.Amjx' 11/11 V? V / G CiR') , 
where \\ ■ stands for the following quantities: 

(i) lfa>0, then \\f\\,/, = {J^, \f\'^^ dxY . 

(ii) If a <Q, let k be the integer part of {—crd) and a = |cr|rf — k be the fractional (positive) part of a. 
Using the standard multi-index notation, where \ri\ ^ rji + . . . + rjd is the length of the multi-index 
V ^ {VIt ■ ■ Vd) G Z'', we define 



ll/lli/.= 



I I \^'^fix)-^'^f{y)\ , „ 
max \d^ f\ = max sup ; ; — IfWc^uw) if a > , 

max sup Id'^fiz)] := H/Hchk") if a = . 

\n\=k zeTSf 



As a special case, we observe that = ||/||cj(R'*)- 

(Hi) By convention, we note ||/||i/o = sup^ggd |/(z)| = ||/||co(R'') = ll/lloo- 

Next we need a regularity result that helps us to compare the C"-norm with the L°°-norm, and we 
combine it with the above interpolation in order to obtain the same rate of decay for all L^-norms. 
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Lemma 5.10 Let m > 0. Let v be the reseated solution to equation (1.3) eorresponding to an initial 
datum uq as in Theorem 2.1 There exists to > 0, a Cz (0, 1) and a constant B > such that v{t, x)~S{x) 
is in C"(r2) and 



(5.22) 



\vit,x)- S(x)\\c^(^n) < B\\v{t,x)- Six)\\^a.^n) V t > io • 



Proof. Let h{t,x) :— v{t,x) — S{x). Since both v and S are solutions to equation (1.3), h solves 
ht^{v- S)t =vt = A(v") +cv = A{{h + S*)™) + c{h + S) . 



By Theorem 2.1 we know that for some to > 0, for any t > to, ||/i(i)||oo can be taken uniformly small 
and V is positive in f2, u = on d^l. The Holder continuity now follows by the nowadays classical results 
of DiBenedetto et al. (cf. the book HJ Chap III, Thm. 1.1 for m > 1 and Chap IV, Thm. 1.1 for 
< m < 1), and holds for a class of equations of the type ht — V ■ A{t, x, h, Vh) + B{t, x, h, V/i), which 
satisfy standard structure conditions: 

A{t,x,h,Wh) ■ Vh > colftr^^lV/ip - M^^t) 
\A{t,x,h,Wh)\ < ci\hr^^\\/h\^ + ipi{x,t) 
\B{t,x,h,Vh)\ < c2\V\hr\^ + ip2ix,t) 

for suitable Ci > and nonnegative ipi. In our case we have that 

Ait^x^h.Vh) ^m{h + Sy^-^Vh and B{t, x, h^Vh) ^ ■ {{h + S)"'-^\' S) + c{h + S) 

clearly satisfy the structure conditions. □ 
The same regularity estimates can be proved for the relative error, at least in the case m > 1. 



5.4 Proof of Theorems 5.6, 5.7 and 5.8 



The result of Theorem 3.5 corresponds exactly to the weighted estimate (5.15) of Theorem 5.6 It 
remains to prove the L°°— estimate (5.19). To this end we combine the results of the previous lemmata. 



h(t,x) := v(t,x) — S{x). Lemma 5.9 gives 



\\h\\oo < A\\hrc4hrp 



1-e 



where we take a the Holder exponent of Lemma 5.10 and we take any p > 0. Then e £ (0, 1) and 
A — A{a,p, oo) > are as in Lemma 5.9 The result of Lemma 5.10 reads 

WHc < B\\h\\^ Vi >to • 
The combination of these two results gives 

(5.23) \\h\\^<A{B\\h\\o.Y\\h\\l-' that is \\v{t,-) - S{-)\\oo < {AB')^ \\v{t,-) - S{-)\\p. 

We now combine the above interpolation inequality with the exponential decay of the weighted L^— norm 
of Proposition (3.5 1: there exists a constant k > such that 

|2 



n 



\v{t,x) - Six)]-" S{x)"'-^ dx 
for alH > ti ^ 1, for all 7 such that 



n 



e{t)\ S^+"'dx < Ke-^(*-*i)£-[^(ti)] 



< 7 < 70 



- 1 



fco(™ 

fci(rn) 



l2 r 



(52Ai)^ \n\ 



2(1 - m) 
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By Holder inequality we have 

\\v{tr)-S{-)\\l 

ll^llll" 



< 



\v(t, x) - S{x)\^ da; < Ke-''^'~''^£[e{ti)] 



so that, combining it with( |5.23 1, for p = 1 we obtain the second inequality |5. 16| of Theorem 5.6 



\\v{t, •) - < (AB^)^ \\v{t, .) - S{-)\\1 = (AB^)^ «;e-^(*-*^)f ^^ic-^^*-*^) 

So far we have concluded the proof of Theorem |5.6| and rescaling back we have proved Theorem |5.7 
It remains to prove Theorem |5.8| that is the PME case m > 1. We just remark that 



\\S\\IZZI-- / (5"") dx- / dist(2;,9f7)-'^dx 
Jn Jn 

the latter quantity being finite for all m > 1, since < {m — l)/m < 1 . □ 



6 Appendix 

6.1 Intrinsic Poincare inequality 

We give a proof of Proposition |3.1[ Notice first that 

(6.1) \2 ( fdx< [ |V/|^da; whenever 

Jn Jn 



= 0. 



We now apply inequality (6.1 ) to the function 

^.g$fda; 



$1 = (5-5*i)*i , 



for which the above orthogonality condition clearly holds. Moreover, we have: 



fdx^ 
n Jn 

[ g^<^ldx 
Jn 



In n 



$1 



dx 



Jn ' In^dx in^ ' 



In addition: 



V/ da; 



In da; 



In^dx 



$f da; 



9^ - 



dx. 



$1 



da; 



= / |Vg$i|' da; 

= [ |Vg$i|' da; 
Jn 



r/o5*?da;1 


7 


ilnnd^l 


Jn 


r/n5*?da;1 


7 


iln^ii^l 


Jn 



Jpg^fda; 
' In^dx 

Jn9^l da; 



' f „ Lg<^idx f 



= / |Vff|'$?da; + Ai / g2^?dx + Ai 

= / |Vg|'$?da; + Ai / g^^^dx-Ai 
Jn Jn 

= / |Vg|2$?da; + Ai / |5-5*j'$?da; 
Jn Jn 



r/nff*?dx1 


7 


Un^?dxJ 


Jn 


7n5*? dx- 


7 


[ /n dx J 


Jn 



/ $^dx-2Ai 
Jn 

dx 



/n*?d^ 



dx 



37 



Summing up we have shown that 



Jn Jn Jn Jn Jn 

which yields the desired inequahty. □ 

We recall some bounds on A2 — Ai. Singer et al. [31], [IS], [3S] proved that for convex domains f7 C M'' 
with diameter diam(f2) and inradius inr(r2), such latter quantity being defined as the supremum of radii 
of balls included in fi: 

diam(f2)2 2 Al < j^^^^^2- 

This bounds can be somewhat improved when further geometrical properties of fl hold, [39] . Notice that, 
by taking O to be a rectangle of sides L and with L large one explicitly computes A2 — Ai — 3n'^/L^, 
and L is close to be the diameter of f2. A lower bound of the form A2 — Ai > Stt"^ / diam{^l)'^ is conjectured 
to be the sharp one. 



6.2 Facts on the Elliptic Problem 



As mentioned above, the stabilization of the solutions w > of the transformed evolution problem ( 1.3 ) 
leads in a natural way to the consideration of the associated stationary solutions, i.e., the solutions of 
the following elliptic problem 

-A(5")=c5 inn, 
S{x) > for X en, 

Six) = for a; € dn, 

where iris < m < 1 and 57 C M.'^ is an open connected domain with sufficiently smooth boundary. Using 
the new variable V — S*™ and putting p = 1/m > 1 the latter problem can be written in the more 
popular semilinear elliptic form 

-AV = cV'P in n, V = on dn 

Note that our restriction m > nis is the exact condition that makes the last problem subcritical, p < Ps- 

• Existence of positive classical solutions 

The question of existence and regularity is well understood in its basic features: 

(a) if < m < 1 for d < 2 or if = < m < 1 for c? > 3, then there exist positive classical solutions 
to equation (see e.g. [7] and references quoted therein, and also [M]). 

(b) if < m < nis and d > 3 then there are cases in which the positive classical solution exists (e.g. if n 
is an annulus) and cases in which it does not exist (e.g., if n is star-shaped) (see e.g. [7] and references 
quoted therein). 

We observe that the geometry of the domain plays a role in the question of existence, but only in the 
subcritical case (b), which is not considered in this paper. Since we assume that m > rrig, the existence 
of at least one positive classical solution is always guaranteed. 

• Uniqueness. In the supercritical case m > TOs considered here, the geometry of n plays a role in the 
uniqueness problem. For example, if d = 1 or if d > 2 and 51 is a ball, then the solution is unique, cf. 
[T]. While when d > 2 and n is an annulus, then the solution is unique only in the class of positive 
radial solutions, cf. [31]. However, there are cases in which the solution is not unique, cf. [3^ I16j . 



Regularity and boundary behaviour. We state now the main bounds for (6.2 1, with explicit constants 



for all 1 < J) < They will give us explicit bounds for the constants ko{p), ki{p) appearing in Theorem 
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4.1 We remark that we already know that ki{p) 1 as p 1, but we have no exphcit bounds for 
them. While providing below such bounds for all 1 < p < Pc, we notice that the resulting estimates will 
not satisfy the above limiting property. The proofs follow by using the arguments that can be found 
for example in [21] for the local bounds, or in |2H1 HO] for the boundary estimates, and they will be 
published separately in [T^]. 

Theorem 6.1 (Local Upper Estimates) Let 51 C M** 6e a bounded domain, and let c > 0. Let V be 

a local weak (sub-) solution in Bji^ C f2 to —AV — cV^ , with l<p<pg=2* — l = {d + '2)/{d — 2). 
Then for any Rao < Rq the following bound holds true: 



(6.2) 



2<j-d(p-l) 

9,- 



for any 



d{p-l) 



< q 



where the constant hi p depends on d, q, p, Ro,Roo o,nd can be explicitly calculated as in |12| . 

Theorem 6.2 (Local Lower Estimates) Let 57 C M"* he a bounded domain, and let c > 0. Let V be 

a local weak solution in Bji^^ d to — AV = cV^ , with I < p < ps = 2* — 1 ^ (d + 2) / (d ~ 2) . Then 
for any e > and for any 



0<q< 



i-2 

2 — 



dLol[eid-l)+s] 



the following bound holds true 
(6.3) 



inf V{x) = ||V||-oo.fl^ > ho^^^^^ 

x6Bh^ iBflol' 



where the constant ho depends on d, q, s, S2, Ro,Roo and can be explicitly calculated as in |12) . 

By means of these upper and lower bounds one can prove quantitative Harnack estimates. 

Theorem 6.3 (Harnack inequality for 1 < p < Pc) Let C M.'^ be a bounded domain, and let c > 
. Let V be a local weak solution in Bj^^ d Q to —Ay — cV^ , with I < p < Pc = d/{d — 2), and assume 
that \\V\\p+i^iig < Kp. Then the following bound holds true for all R^o < Ro- 



(6.4) 



sup V{x) < Hp inf V{x) 



where the constant "Hp depends on d, p, S2, Rq, Rooi I^p o-nd can be explicitly calculated as in [12' . 

We now compare solutions corresponding to different p and c, and this can be done for 1 < p < Pc , 
since we need the quantitative Harnack inequalities of Theorem |6.3[ that hold only in that range of p. 



We recall that we are now choosing c = 1/[(1 — m)T] = p/[{p— 1)T] so that by Proposition 5.4 we have 
that c Ai as p — > 1. Hence Hp in the above Theorem has a finite limit Hi as p 1. 



Theorem 6.4 Let Up be a weak solution to the Elliptic problem 4.1 without any assumption on Xp > 0, 
and let 1 < p < Pc = d/{d — 2). There exists a 6{Q) = 6 > independent of Up, such that 



(6.5) 

where 

(6.6) fc„(p) 



1 hi^p 
2'Hp hi I 



Ri 



d-2 



klip) :=2Hi 



hi.p 
hi.i 



Ri+S 



-id-2) 



and Cp are the constants in the upper bounds of Theorem 6.1 and Hp is the constant in the Harnack 
inequality of Theorem\6.3\ 



39 



6.3 Maximum and comparison principles on small sets 

The maximum and comparison principle do not hold in general for solutions to nonlinear elliptic equa- 
tions. This is an important characteristic of elliptic equation in general and does not necessarily depend 
on the nonlinearity. Indeed in the linear case, if one consider the Dirichlet problem for the equation 
—Alt = Xu with A > All it happens for instance that for A = A2 > Ai > the corresponding second 
eigenfunction $2 has at least a change of sign, hence no global maximum nor comparison principle is 
allowed to hold. 

In any case, we can still prove a (local) maximum and comparison principle on small sets: we are 
going to extend to our framework an idea originally due to Serrin, see for example the book |37j where 
this idea is applied here to a different class of nonlinear elliptic equations. We just state the Theorem 
here, a complete proof will appear separately in [12j . 

Theorem 6.5 (Comparison with supersolutions on small sets) Let B C Mf'' be a bounded con- 
nected domain, let p > 1, X > and 

{-Au = XuP inB 
-Au>XuP inB 
u> u on dB 

0<u,u< M inB 

and assume that \B\ < uj^j (2p X M^^^^'^ . Then, we have that u > u in B. 
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